Orbital Ordering in Paramagnetic LaMnOs and KCuF 3 . 
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Ab-initio studies of the stability of orbital ordering, its coupling to magnetic structure and its 
possible origins (electron-phonon and/or electron-electron interactions) are reported for two per- 
ovskite systems, LaMnOs and KCUF3. We present a new average spin state calculational scheme 
that allowed us to treat a paramagnetic state and to succesfully describe the experimental mag- 
netic/orbital phase diagram of both LaMnOs an d KCUF3 in crystal structures when the Jahn- Teller 
distortions are neglected. Hence, we conclude that the orbital ordering in both compounds is purely 
electronic in origin. 
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It is known from earlier studiesEro that magnetic order- 
ing is coupled to the orbital ordering (00) - magnetic 
interactions depend on the type of occupied orbitals of 
transition-metal ions. (A very detailed overview on mag- 
netic and orbital ordering (00) in cuprites and man- 
ganites was given recently by Oles et aB.) For exam- 
ple, the 00 in LaMn03 reduces the ferromagnetic (FM) 
contribution from the e g orbitals to the interlayer ex- 
change coupling and is responsible for the stability of 
the magnetic ground state structure (which is A-type 
antiferromagnetic (AFM) below Tjv=141 Ka) — as de- 
scribed both|-|experimentally and theoretically on this 
manganiteo'QU. In_KCuF3, the 00 was obtained in a 
spin-orbital modelo as a result of an exchange interac- 
tion which correctly describes both orbital and spin align- 
ments simultaneously; thus the spin and orbital degrees 
of freedom are mutually coupled. 

A more complicated picture of the fundamental con- 
nection between 00 phenomena and forms of the spin 
alignment follows from experimental measurements that 
characterize an 00 parameter and magnetic interactions 
in transition-metal oxides. Using the diople resonant x- 
ray scattering technique, Murakami et alU found a sharp 
disappearence of 00 at much higher temperatures (~780 
K) than Tjv in LaMn03. They suggested a coupling of 
the spin and orbital degrees of freedom due to the small 
decrease of the order parameter of the orbital structure 
above Tjy, even despite the presence of 00 in the para- 
magnetic (PM) phase of LaMn03. The most-striking re- 
sult was obtained for bilayered LaSr 2 Mn207l3, namely, a 
competition between the A-typ£ AFM spin ordering with 
Tat «170 K and the CE-typeo charge/orbital ordering 
(COO), which exists between T/v and COO transition 
temperature, Tcoo=210 K. With decreasing tempera- 
ture, the development of the COO phase is disrupted by 
the onset of the A-type AFM ordering at Tjv- 

In this paper, we investigate the stability of observed 
OO's without appealing to magnetic interactions by mod- 
eling the PM phase of two manganites, LaMnOs and 
KCUF3. To this end, we examined all magnetic con- 
figurations (MC's) that are mathematically possible in 
a unit cell of the perovskites which consists of 4 Mn 
(or Cu) atoms. There are only eight such spin sets - 
as schematically presented in Fig. 2(a). They include 



the possibility of having FM (labelled "1"), and C-type 
("3"), A-type ("6") and G-type ("7") AFM spin order- 
ings - as well as four MC's that have a mathematical 
meaning, but cannot be realized as an individual mag- 
netically ordered solution. Then, to model a PM state, 
during the iterations towards self-consistency, we average 
the orbital occupation matrices (that contain information 
on the orbital polarization and 00 in the unit cell) and 
potential parameters of all atoms over eight spin align- 
ment configurations after each iteration. As a result, we 
have an averaged spin state and can analyze the stabil- 
ity/presence of orbital polarization in such a modeled 
PM phase. Self-consistency of the averaged spin state 
calculation (ASSC) is determined by convergence of the 
averaged total energy, charge densities and orbital occu- 
pation matrices. 

Typically, the origin of the orbital polarization is be- 
lieved to be the electron-phonon interaction (the Jahn- 
Teller (JT) distortion) and the electron-electron interac- 
tion. A cooperative JT effect which sets in well above 
the magnetic transition temperature, stabilizes the par- 
ticular order of the partly filled e g orbitals that arg 
close to orbital degeneracya. Monte Carlo investigations^ 
showed a stabilization of the experimental magnetic and 
00 in LaMnOs only after including lattice distortions. 
However, a spin-orbital model for insulating undoped 
LaMn03 derived inO qualitatively explains the observed 
A-type AFM ordering as stabilized by a purely electronic 
mechanism. 

We investigated two typical pseudocubic perovskites: 
(i) LaMnOs, in which the strong JT effect breaks the de- 
generacy of the electronic configuration of Mn 3+ (t^g^g) 
and directly affects the e g orbital population. A recent 
neutroiwaowder diffraction study, together with thermal 
analysial3, showed that LaMn03 undergoes a structural 
phase transition at Tjt=750 K above which the 00 dis- 
appears. To establish the origin of 00 in LaMnOs, we 
performed ASSC's for the high temperature crystal struc- 
ture-wjacrc JT distortions are absent. Since diffraction 
dataLLj'til can be interpreted using two crystal symme- 
tries - double-cubic orthorhombic Pbnm structure and 
rhombohcdral i?3c crystal structure - we did calcula- 
tions for both, (it) Another example of a perovskite com- 
pound with partly occupied e g orbitals resulting in strong 
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JT distortions and strong one-dimensional AFM interac- 
tions iSpKCuF 3 . There are two tetragonal crystal types of 
KCuFalla, which are stable in a wide temperature range, 
and their A-type AFM ordering,temperatures are equal 
to 20 K and 38 K. It was shownEj for this compound that 
the ground state with the lattice distortion and the or- 
bital polarization (and ordering) can be stabilized only by 
taking into account the Coulomb correlation between d- 
shell electrons (LDA+U approaclJiil) . We use Coulomb 
and exchange parameters, U=8 eV and J=0.88 eV for 
both compotHids as those relevant to the 3d transition 
metal oxideai-l'E-a 

LaMnOz- As the first step, we performed LDA+U 
calculations for the experimental A-type AFM phase of 
LaMnC"3 with low-T Pbnm structure (where strong JT 
distortions are present). The self-consistent non-diagonal 
occupation matrix for the spin density of the e g subshell 
(3z 2 — r 2 and x 2 — y 2 orbitals) of one particular Mn atom 
is found to be 

'W 'W ^ o.29 0.59 J ' 

Its diagonalization gives two new e g orbitals: 4>i — 
3x 2 — r 2 (3y 2 — r 2 for the second type of Mn atom) with an 
occupancy of 0.81 and = z 2 ~y 2 (z 2 —x 2 ) with an occu- 
pancy of 0.20. Since a Mn 3+ ion has formally one electron 
in the partially filled subshell and the t 2g subshell is to- 
tally filled, we presented the resulting OO by plotting the 
angular distribution of the e g spin density in Fig. 1(a): 
P(M) = E m „>W -n-L n ,)Y m (6,4)Y m ,(8,cb), where 
the Y m (9,<p) are corresponding spherical harmonics. 

The resulting picture of the OO in A-type AFM 
LaMnOa is in very good agreement with that given by 
GoodenoughEl and recently detected using resonant x-ray 
scattering!! Based on this result, we can now move to 
the case of T > TV to investigate the region of the mag- 
netic/orbital phase diagram where the OO exists in the 
PM phaseQ. To this end, we performed the model ASSC's 
described above for the high-T (at T=798 K) structure 
of LaMnOa with experimental lattice parameterst2l. The 
PM iterations were started from a uniform distribution 
of e g electrons over the 3z 2 — r 2 and x 2 — y 2 orbitals of 
all Mn atoms. In Fig. 1(b), we plot the "total energy" of 
the PM state of the high-T structure of LaMn03 (crosses) 
obtained by averaging the total energies of the eight pos- 
sible MC's. This total energy curve becomes saturated 
after the ~20th iteration; hence, we may call this ASSC 
result self-consistent. 

To illustrate the stability of OO in the PM phase and 
to compare the orbital polarization with that calculated 
for the experimental A-type AFM phase, we present the 
OO as a difference between the occupancies of the ("diag- 
onalized") <\>\ and cj>2 orbitals. In Fig. 1, this difference 
is shown plotted against the number of PM iterations 
(crosses). We found that starting with equal 2>z 2 — r 2 
and x 2 — y 2 orbital occupancies (and with a diagonal 
occupation matrix), the value of the non-diagonal ele- 



ments of the occupation matrix grows during the ASSC 
leading to a saturated orbital polarization of the right 
type (the same obtained for A-type AFM LaMn03, Fig. 
1(a)) around the tenth iteration (Fig. 1). Note that 
the OO in the PM state is the only stable solution. Forc- 
ing different orbital polarizations (ferro-orbitally ordered, 
etc.) as starting ones for the PM calculation resulted 
in the type of OO obtained above by changing the cor- 
responding (non-diagonal) orbital occupancy during the 
iterations towards self-consistency. The horizontal line in 
Fig. 1 stands for the difference between cf> 1 and <j>2 orbital 
occupancies obtained from the usual spin-polarized self- 
consistent LDA+U calculation for the A-type AFM phase 
of LaMnOs in the low-T Pbnm structure. Comparing the 
orbital polarizations of the PM and AFM phases (Fig. 1), 
we conclude that the AFM-PM transition in LaMn0 3 re- 
sults in a slight decrease of the difference between the <j>\ 
and <p2 orbital occupancies but does not supress the OO. 
This finding agrees with the observed behavior of the OO 
parameter when T increases (T > T^p. 

We need to point out here one essential feature of the 
ASSC. Since we consider a spin state averaged over eight 
MC's and therefore an average ocupation matrix of such 
a PM state, it is interesting to analyze the deviation of 
the orbital polarization in each of these MC's from the 
average value. As seen from Fig. 3(b), where we plot 
the and <p2 orbital occupancy difference for each of 
the eight MC's after the first (broken line) and fiftieth 
(solid line) PM iteration, the deviations from the aver- 
age value become smaller with iterations - they decrease 
from 36% to 12%, correspondingly. (Further increases 
of the iteration number do not supress the deviations.) 
The variations at the beginning of the ASSC and for the 
result obtained, c.f., Fig. 2(b), have a close analogy in 
oscillation behavior, i.e., the location of the e g orbital 
occupancy difference stays higher (lower) with respect to 
the average value for the same particular MC during the 
iterations. Hence, one can separate the eight MC's into 
three groups: the lowest (MC's 1 and 2), the highest (4 
and 7) and those close to the average value (3,5,6 and 8). 

We analyzed possible reasons for this separation by cal- 
culating the number of FM or AFM neighbours aligned 
to the first Mn (whose spin does not change in all eight 
MC's, Fig. 2(a)). First, we found that the behavior of 
the e g orbital occupancy difference curves (Fig. 2(b)) is 
determined by the nearest neighbours' spin alignment. 
For both the 1st and 2nd MC's, all six nearest neigh- 
bours of the first Mn have the same spin direction, so 
FM couplings with all nearest neighbours give the small- 
est <j>x — cf>2 value. By contrast, the largest contributions 
were obtained from MC's 4 and 7, which have six AFM 
couplings between the first Mn and its nearest neigh- 
bors. The third group of MC's possesses two FM and 
four AFM couplings (MC 3 and 8), giving almost an ex- 
act average value <\>\ — <j>2, and four FM and two AFM 
couplings (MC 5 and 6), both lying below the average 
value. (Note here, that it would be sufficient to consider 
only four MC's: the 1st (which corresponds to the FM 
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phase), the 3rd (C-type AFM), the 6th (A-type AFM) 
and the 7th (G-type AFM). Averaging over these four 
MC's underestimates the value averaged over all eight 
MC's by only 0.7 and 0.1 % at the first and 50th PM 
iteration, respectively.) Thus, we can conclude that the 
AFM spin alignment seems to be more preferable for 00 
than the FM one; however, the decrease of the </>i — <f>% 
deviation against the average values (from 47 to 17 % 
for the first and 50th iteration, respectively) in the FM 
phase also supports the importance of the FM coupling 
in ASSC. 

Let us now discuss some peculiarities in the crystal 
structure of LaMnOs related to the JT distortioaSj. As 
found by recent neutron-powder diffraction studies^, the 
JT transition in this perovskite compound occurs at 
T=750 K. If we define the degree of tetragonal distor- 
tion as Sjt — (di — d s )/(di + d s )/2, where di and d s 
denote the long and short Mn-0.-bond distances, then 
according to the structural datal!3, Sjt changes with 
heating from <5jt=0.133 for the room-T Pbmn struc- 
ture to (5jt=0.023 for the crystal structure at 798 K 
(which is also orthorhorobjc Pbmn, but described as 
"double-cubic" perovskiteta) . Based on the results ob- 
tained above, we may conclude that strong lattice distor- 
tions do not influence the OO. However, Sjt is not zero in 
the high-T structure and we have slightly different Mn-0 
bond distances which break the cubic symmetry, and so 
the JT effect may still play some role in 00. To check this 
possibility, we performed ASSC for the rhombahjsdral 
i?3c structure with lattice parameters taken fron£3. All 
Mn-0 distances are equal in this structure, so <5jt=0.000. 
For this case, the calculated difference between occupan- 
cies of the <f>i and (j>2 orbitals against the number of PM 
iterations is shown in Fig. 1 (circles). This difference 
is smaller than those obtained for the two orthorhombic 
structures; nevertheless, we conclude that the 00 does 
not disappear in regular Mn06 octahedra, and hence JT- 
distortions are not the origin of the 00 in LaMn03. 

KCuF^. As shown by Kugel and KhomskiiB, KCUF3 
is an example of a system in which the exchange interac- 
tion alone results in the correct 00. Ab-initio LDA+U 
investigations^ confirmed the electronic origin of the or- 
dering: the coupling to the lattice is not a driving force 
for the orbital (and magnetic) ordering, but the lattice 
follows the orbital state. Hence, we performed ASSC's 
for a model structure of KCuF 3 in which cooperative JT 
lattice distortions were neglected. We used the tetragonal 
(PA/mmm) crystal structure with a=5.855 A, c=7.846 
A and the coordinates: K (0, 0.5, 0.45), Cu (0, 0, 0), Fl 
(0, 0, 0.45) and F2 (0.25, 0.25, 0). In this structure, the 
CuFg octahedra are slightly compressed along the c axis: 
the distance from the Cu atom to the apical Fl atom 
D(Cu-Fl) = 1.96 A, while in the afe-plane D(Cu-F2)=2.07 
A (without the quadrupolar deformation in the ab plane 
that is present in the experimental structure). 

Since this perovskite compound possesses two ener- 
getically equivalent types of OO's characterized pb,y al- 



performed four ASSC's which started with different e g - 
electron distributions over Cu e g orbitals: the orbital po- 
larizations chosen corresponded to F-type ferro-orbital 
and A-, C- or G-type antiferro-orbital alternations of 
x 2 — z 2 /y 2 — z 2 orbitals (shown in the upper row of OO's 
in Fig. 3); the C- andpG- orbital configurations were 
observed experimentallyta. Since a Cu 2+ ion (d 9 con- 
figuration) has one hole in the e g state, we again rep- 
resent the OO's by a 2x2 e 9 -orbital occupation matrix, 



As a start, we used 



0.53 
±0.28 



whose diagonalization gives x 2 



±0.28 
0.15 



2 ) orbitals. 



In contrast to LaMnOs, we found some stable solutions 
with different OO's in KCuF 3 . In Fig. 3, we present 
the total energy of the PM phase for each of the orbital 
ordered states obtained against the number of PM iter- 
ations. The corresponding angular distributions of the 
e g spin density for every solution (at the 200th iteration) 
are shown (bottom row of OO's in Fig. 3); the antiferro- 
orbitally ordered state (A) has the highest total energy. 
Although it seemed to be stabilized with the lowest total 
energy and with almost fully occupied 3z 2 — r 2 orbitals at 
all Cu sites at around the 70th PM iteration, further PM 
iterations resulted in a sharp increase of the total energy 
of this A-type solution, and at the 200th iteration the OO 
is an alternation of 3z 2 — r 2 and 3x 2 — r 2 /3y 2 — r 2 orbitals. 
For the calculation, which was started from F-type ferro- 
orbital ordering, the same x 2 — y 2 orbitals at all Cu sites 
(upper row of Fig. 3) change to almost fully occupied 



3z 2 



orbitals during the PM iterations (bottom row 



in Fig. 3). 

More preferable in energy are two (almost degenerate) 
solutions which correspond to G- and C-type OO. Their 
4>i and <p2 orbitals, x 2 — z 2 and y 2 — z 2 , stagger in the 
ab plane, but have ferro- or antiferro-orbital alignments 
along the c axis for C or G-type ordering, respectively. 



Btical and experi- 
Thus. from the 



This result agrees with previous . 
mental studies on OO in KCUF3L 
behavior of the total energy curves of these four solu- 
tions, we conclude that PM KCUF3 with uniform CuFg 
octahedra possesses two OO's of the right type. 

In summary, the new ASSC scheme presented to treat 
a PM state succesfully described the experimental mag- 
netic/orbital phase diagram of both LaMnOs and KCUF3 
without inclusion of JT distortions. The OO in the 
PM phase is found to be the same type as in the AFM 
phase. A small decrease of the order parameter of the 
OO upon the AFM-PM transition is in agreement with 
the observed one. Finally, we conclude that OO in both 
LaMnOs and KCUF3 is of purely electronic origin. 

Work supported by the Russian Foundation for Basic 
Research (grant RFFI-01-02-17063) and the U.S. Depart- 
ment of Energy (grant No. DE-F602-88ER45372). 



tcrnation of the 



and y 



z 1 orbitalsB'lHi we 



3 



1 J. B. Goodenough, Phys. Rev. 100, 564 (1955). 

2 K. I. Kugel and D. I. Khomskii, Sov. Phys. Usp. 25, 231 
(1982); Zh. Eksp. Teor. Fiz. (Sov. Phys. JETP) 64 (37), 
1429 (1973). 

3 A. M. Oles, M. Cuoco, and N. B. Perkins, in: Lectures on 
the Physics of Highly Correlated Electron Systems IV, AIP 
Conference Proceedings, edited by F. Mancini (New York, 
2000), Vol. 527, p. 226-380. 

4 E. O. Wollan and W. C. Koehler, Phys. Rev. 100, 545 
(1955). 

5 K. Terakura et al, Mater. Sci. Engineering B63, 11 (1999). 

6 Y. Murakami et al, Phys. Rev. Lett. 81, 582 (1998). 

7 T. Chatterji et al, Phys. Rev. B 61, 570 (2000). 

8 T. Hotta, S. Yunoki, M. Mayr and E. Dagotto, Phys. Rev. 
60, R15009 (1999). 

9 L. F. Feiner and A. M. Oles, Phys. Rev. 59, 3295 (1999). 

10 J. Rodriguez-Carvajal et al, Phys. Rev. B 57, R3189 
(1998). 

11 Q. Huang et al, Phys. Rev. B 55, 14987 (1997). 

12 M. T. Hutchings, E. J. Samuelsen, G. Shirane and K. Hi- 
rakawa, Phys. Rev. 188, 919 (1969); A. Okazaki, J. Phys. 
Soc. Jpn. 26, 870 (1969); 27, 518 (1969). 

13 A. I. Liechtenstein, J. Zaanen and V. I. Anisimov, Phys. 
Rev. B 52, R5467 (1995). 

14 V. I. Anisimov, J. Zaanen and O. K. Andersen, Phys. Rev. 
B 44, 943 (1991); V. I. Anisimov, F. Aryasetiawan and 
A. I. Liechtenstein, J. Phys.: Cond. Matter 9, 767 (1997). 

15 S. Satpathy, Z. S. Popovic and F. R. Vukajlovic, Phys. 
Rev. Lett. 76, 960 (1996); I. Solovyev, N. Hamada and K. 
Terakura, Phys. Rev. B 53, 7158 (1996). 

16 A. Vegas et al, Acta Crystallogr., Sect. B: Struct. Sci. 42, 



167 (1986). 

17 J. F. Mitchell et al, Phys. Rev. B. 54, 6172 (1996). 

18 K. Hirakawa and Y. Kurogi, Progr. Theor. Phys. (Suppl.) 
46, 147 (1970). 



FIG. 1. Dependence of the <f>\ — <j>2 difference on the num- 
ber of PM iterations for the high-T PM phase, <5j T =0.023 
(crosses) and R3c PM phase, <5 JT =0.000 (circles) of LaMn0 3 . 
The horizontal line stands for calculated OO of Pbmn A-type 
AFM phase, Sjt~0.133. In the insert, (a) denotes the angu- 
lar distribution of the e 9 -electron spin density in the Pbmn 
A-type AFM phase of LaMnOs from the LDA+U calculation 
and (b) the total energy averaged over all MC's as a function 
of the number of PM iterations for the two PM phases. 

FIG. 2. (a) Schematic illustration of all possible MC's in 
the unit cell which contains four formula units; only the ac 
plane is shown, (b) Deviations of the cj>i — cj)2 difference in each 
of the MC's considered from the averaged value as a result of 
ASSC's of the high-T PM phase of LaMn0 3 . Only the first 
(broken line) and fiftieth (solid line) iterations are shown. 

FIG. 3. Total energies of four solutions with different OO's 
for KCUF3 as a function of the number of PM iterations. The 
corresponding OO (F-type ferro-orbital and A-, C- or G-type 
antiferro-orbital orderings) are shown for the zeroth (top row) 
and the 200th iteration (bottom row). 
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